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Abstract 


Double-barrier  resonant  tunneling  structures  often  contain  undoped  spacer 
layers  separating  heavily  n -doped  regions  from  active  regions  of  the  structure. 
Conduction  band  electrons  coming  from  the  heavily  doped  regions  diffuse  into  the 
undoped  regions,  producing  potential  energy  profiles  that  can  confine  hole  states 
within  the  spacer  layers.  A  classical  model  is  used  to  predict  the  potential  energy 
profiles  in  these  structures,  and  the  associated  hole  binding  energies  and  wave 
functions  are  calculated  numerically.  (  ) 
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1.  Introduction 


Double-bamer  resonant  tunneling  (DBRT)  structures  (and  in  fact  a  variety  of 
different  semiconductor  heterostructures)  often  contain  heavily  doped  re¬ 
gions  and  undoped  “active”  regions  separated  by  an  undoped  spacer  layer. 
Attempts  to  model  the  electronic  states  in  these  structures  always  involve 
assumptions  regarding  the  distribution  of  extended-state  electrons  (or  holes) 
near  the  abrupt  doping  profile  edge.1  For  example,  n-type  DBRT  structures 
often  contain  undoped  spacer  layers  (about  100  A  thick)  surrounding  the 
bamer/quantum-well  (QW)  region  of  the  structure.  These  layers  are  in  turn 
surrounded  by  heavily  n-doped  regions  that  facilitate  electrical  contact  to  the 
sample.  The  present  paper  deals  with  the  nature  of  the  band-bending  near 
spacer  layers  and  the  hole  states  that  are  confined  (in  one  spatial  dimension) 
by  this  space-charge  potential. 

We  use  a  classical  model  of  the  charge  distribution  to  derive  expressions  for 
the  potential  profile  in  unbiased  DBRT  structures.  We  then  integrate  a  one- 
band  effective  mass  Schrodinger  equation  through  the  region  of  interest  to 
obtain  the  hole  states  confined  by  the  potential.  In  a  concluding  section,  we 
comment  on  the  classical  model’s  weakness  and  indicate  how  the  theory 
might  be  improved  so  that  the  potential  profile  and  binding  energy  estimates 
are  more  reliable. 

2.  Model  and  Calculations 

Consider  the  energy  band  diagram  of  a  typical  DBRT  structure,  shown  in 
figure  1.  The  host  semiconductor  might  be  GaAs  and  the  barrier  layers 
AlGaAs.  We  assume  that  the  region  between  the  vertical  dotted  lines  is 
undoped  while  the  region  outside  the  dotted  lines  is  uniformly  doped.  We 
further  assume  for  simplicity  that  all  donor  atoms  are  ionized  and  that  the 
electrons  associated  with  them  occupy  extended  conduction  band  states.  In 
modeling  the  free  electron  distribution  in  this  structure,  we  replace  the 
discrete,  randomly  distributed  ionized  donor  atoms  by  a  uniform  positive 
background,  and  we  ignore  potential  fluctuations.  We  now  ask,  how  do  the 
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Figure  1.  Energy  band 
diagram  of  a  double¬ 
barrier  resonant 
tunneling  structure. 
Region  between  vertical 
dotted  lines  is  undoped. 
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electrons  distribute  themselves  in  the  presence  of  this  background  charge 
density?  We  address  this  question  from  a  classical  point  of  view,  commenting 
later  on  improvements  to  the  model. 

At  T  =  0  K,  the  electron  charge  density  will  perfectly  match  the  positive 
background  charge  density,  and  the  electrostatic  potential  will  be  zero 
everywhere.  At  a  finite  temperature,  some  electrons  will  diffuse  into  the 
undoped  region  (on  a  distance  scale  comparable  to  a  Debye  length)  leading 
to  dipole  layers  centered  near  the  doping  edges.  This  will  produce  an 
electrostatic  potential  that  decreases  as  one  enters  the  undoped  region  from 
either  side  of  the  structure  (fig.  1  corresponds  to  this  classical  picture  at  zero 
temperature).  Let  z  =  0  correspond  to  the  doping  edge  on  the  left  side  of  figure 
1  (vertical  dotted  line),  with  z  increasing  towards  the  double-barrier  structure. 
We  take  the  spacer  layer  to  have  a  thickness  t  and  assume  that  electron  density 
is  zero  inside  the  barrier  and  QW  regions  of  the  structure.*  This  implies  that 
electrons  which  diffuse  past  the  doping  edge  must  reside  in  the  spacer-layer 
region,  and  therefore,  the  electric  field  at  z  =  t  is  zero. 

In  the  classical  approach  to  the  problem,  the  electron  number  density  n{z)  is 
given  by 

n(z)  =  nd  exp[/fe0(z)]  ,  (1) 

where  0=  1/kT,  nd  is  the  donor  density,  e  is  the  electronic  charge,  and  <p(z)  is 
the  electrostatic  potential  (assumed  equal  to  zero  at  z  =  -~).  Expression  (1) 
is  readily  obtained  by  assuming  that  (a)  0(z)  changes  slowly  enough  that,  at 
any  point  z,  the  density  of  conduction  band  states  is  what  one  would  have  in 
a  bulk  semiconductor  (without  an  electric  field)  whose  band-edge  has  been 
shifted  by  the  local  value  of  the  potential  energy  (-e^(z)),  (b)  the  chemical 
potential  satisfies  the  inequality  Ec  -  e<p(z)  -  fi  »  kT,  where  Ec  is  the 
conduction-band-edge  energy,  and  (c)  n(z)  is  always  much  greater  than  the 
intrinsic  charge  density  in  the  conduction  band.  The  first  two  assumptions 
will  not  be  satisfied  for  some  choices  of  nd,  T,  and  t.  In  particular,  assumption 
(b)  is  certainly  violated  over  some  region  of  the  structure  for  degenerately 
doped  materials.  Nevertheless,  corrections  due  to  election  degeneracy  are 
sometimes  minor  and,  as  a  first  estimate  of  the  potential  profile  and  hole¬ 
binding  energies,  the  nondegenerate  result  will  suffice.  Next  we  impose  the 
requirement  that  Poisson’s  equation  be  satisfied: 

d2<t>(z)  4  Jtp(z)  (2) 

dz2  e  ’ 

where  e  is  the  static  dielectric  constant  of  the  host  semiconductor,  and  p(z)  is 
given  by 

e[rid -  n(z)]  fbrzcO 

p(z)  = 

-en(z)  forz>0 

*We  assume  that  the  barrier  height  is  infinite. 
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By  combining  equations  (1),  (2),  and  (3),  one  obtains  the  equations 

fbrz>0,  and 


d  V(z)  \  r  rrr  X1 
- —  =  -  ^  exp[-(/(z)] 


dz * 


ll 


(4) 


d2Ujz) 
dz 2 


=+jL(l-expK/(z)]) 


forz<0. 


(5) 


where  U(z)  is  a  dimensionless  potential  energy  defined  by  U(z)  =  -/3 e<P(z),  and 
/  is  the  classical  Debye  length  given  by 


/ . 

\4  jcnde2) 


(6) 


Equation  (4)  can  be  integrated  analytically  to  obtain,  for  the  z  >  0  region. 


U{z)  =  U(t)  +  2  In 


cos  0o 


exp[-E/(fl/2]z ' 

ill 


(7) 


where  0O  is  defined  by  cos20o  =  exp([/(0)  -  U(t)). 

The  solution  to  equation  (5)  can  be  given  in  terms  of  a  definite  integral  (which 
cannot  itself  be  profitably  expressed  in  terms  of  special  functions).  We  obtain, 
for  z  <  0, 

(/(z)  =  G-1jM[(/(0)]+^)  .  (g) 

The  single-valued  function  G  is  given  by 


G(x)  =  Af(x)  +  i2  ln[x/(/(0)]  , 


where  Mix)  is  defined  by 


M{x)  = 


dy 


1 

.Vy  +  exp(-y)  -  1 


y 


(9) 


The  solutions  obtained  for  differential  equations  (4)  and  (5)  both  require  the 
specification  of  boundary  conditions.  We  have  already  implicitly  imposed 
continuity  of  U(z)  at  z  =  0.  By  requiring  that  the  electric  field  at  z  =  0  (which 
is  proportional  to  U'iO))  be  continuous,  we  obtain,  equating  first  integrals  of 
(4)  and  (5), 


exp[-C/(r)]  =  1  - 1/(0)  . 


(10) 


Another  relation  between  C/(0)  and  Uit)  is  obtained  by  evaluating  equation  (7) 
at  z  =  t.  After  combining  this  second  relation  with  equation  (10),  we  obtain  a 
transcendental  equation  for  Uit),  namely. 


1  -  Uit)  -  exp  [-C/(01  =  2  Inf cosj^ij^lj 

\  i2l 


(11) 
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This  expression  contains  only  one  unknown,  U(t),  which  can  be  found  by 
root-solving.  Once  U(t)  is  known,  U( 0)  can  be  obtained  from  (10).  The  results 
can  then  be  used  in  the  solutions  to  (4)  and  (5)  to  calculate  the  potential  energy 
function. 

We  have  carried  out  this  procedure  and  in  figure  2a,  we  plot  the  valence-band- 
edge  diagram  of  a  DBRT  structure  at  three  temperatures:  T  =  77, 150,  and 
300  K.  Figure  2b  shows  the  depth  of  the  hole-confining  potential  well  as  a 
function  of  temperature.  We  have  chosen  the  structure’s  parameters  as  fol¬ 
lows:  t  =  150  A,  nd  =  2  x  101S  cm-3, 150-meV  AlGaAs  barrier  height,  42.25-A 
QW  width,  and  62.5-A  barrier  width  (the  latter  three  are  irrelevant  for  the 
calculation  of  the  potential).  We  now  integrate  a  one-band  effective  mass 
Schrodinger  equation  through  the  spacer-layer  region  of  the  structure.  In 
figures  3a  and  3b,  we  plot  the  T=  300  K  hole-confining  potential  (upside  down 
and  reversed)  along  with  the  confined  heavy-  and  light-hole  wave  functions, 
respectively.  Figures  3c  and  3d  contain  similar  results  at  T  =  77  K.  The  wave 
functions  are  labelled  by  the  binding  energies.  These  results  were  obtained  by 
numerically  integrating  a  one-band  effective-mass  Schrodinger  equation  for 
i//(jc)  (jc  =  t  -  z),  subject  to  the  requirement  that  y{x)  go  to  zero  far  inside  the 
barrier  region  (x  «  0)  and  far  inside  the  doped  region  (x  » t),  and  that  y/Vm* 
be  continuous  at  the  barrier  interface  (where  m*  is  the  hole  effective  mass). 
In  our  calculations,  we  used  =  0.34  m0  and  =  0.094  m0  for  the  heavy- 
and  light-hole  masses  in  the  host  semiconductor,  and  =  0.433  m0  and 
mtt=  0.109  m0  for  heavy-  and  light-hole  masses  in  the  barrier  layer  (values 
which  correspond  to  an  A1  fraction  of  about  30  percent). 

We  believe  that  if  electron-hole  pairs  were  generated  within  a  diffusion  length 
of  the  doping  edge,  holes  could  readily  migrate  into  the  spacer  layer  region, 
occupying  the  lowest  hole  subband  calculated  above.  One  would  then  expect 
to  see  luminescence  radiation  resulting  from  transitions  between  extended- 
state  conduction  band  states  and  the  lowest  hole  subband,  and  (maybe) 
between  conduction  band  DBRT  resonance  states  and  the  lowest  hole 


(a)  (b) 


Figure  2.  (a)  Valence-band-edge  energy  diagram  of  double-barrier  resonant  tunneling  structure  at  T  =  300, 
150,  and  77  K.  Spacer  layers  are  150  A  thick,  and  doping  density  of  ft -doped  regions  is  2  x  10u  cm4,  (b)  Hole- 
confining  potential  well  depth  versus  temperature  for  structure  with  same  parameters. 
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Figure  3.  Hole-confining  potentials  (upside  down)  plotted  with  bound  heavy-  and  light-hole  wave  functions: 
(a)  heavy-  and  (b)  light-hole  states  at  T  =  300  K;  (c)  heavy-  and  (d)  light-hole  states  at  T  =  77  K. 


subband.  However,  seeing  such  features  might  prove  difficult.  The  larger 
binding  energies  (which  are  most  easily  separable  from  direct-gap-related 
features  in  luminescence)  occur  in  the  deeper  wells,  which  in  turn  occur  at  the 
higher  temperatures  (according  to  the  classical  picture  described  above).  On 
the  other  hand,  the  higher  the  temperature,  the  broader  are  the  photolumines¬ 
cence  line  shapes,  particularly  in  heavily  doped  samples.  We  believe  that  one 
should  see,  at  least,  a  luminescence  band  associated  with  transitions  between 
the  low-energy  conduction  band-edge  states  and  the  confined  lowest-sub- 
band  heavy-hole  state.  We  have  not  calculated  the  expected  position  (or 
shape)  of  such  a  band  as  a  function  of  temperature  but  believe  that  it  mignt 
have  a  characteristic  temperature  dependence  arising  from  the  temperature 
dependence  of  the  quantum  confining  potential  in  the  spacer  layer.  We  have 
made  some  preliminary  piezoreflectance  measurements  on  a  DBRT  sample 
(with  parameters  identical  to  those  used  in  the  model),2  but  have  not  seen  any 


1  R.L.  Tober  and  J.  D.  Bruno,  Proc.  SPIE 1286  (in  press). 
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convincing  evidence  of  the  quantum-confined  hole  states  in  the  spectra.  The 
experimental  effort  is  still  in  progress. 

3.  Conclusions 

The  assumptions  upon  which  the  classical  model  relies  are  not  well  satisfied. 
We  have  ignored  quantum-mechanical  effects  in  determining  the  distribution 
of  conduction-band  electrons.  We  then  used  the  very  potential  obtained  from 
the  resulting  electron  distribution  in  an  effective-mass  Sc hrodinger  equation 
to  determine  quantum-confined  hole  states.  Our  rationale  fordoing  this  rested 
on  the  lack  of  a  confining  potential  for  electrons.  We  therefore  assumed  that 
the  local  density  of  electron  states  was  modified  less  by  the  potential  than  was 
the  local  density  of  hole  states.  This  rationale,  however,  is  not  very  satisfying. 
Furthermore,  one  can  expect  that  even  at  zero  temperature,  where  the  classical 
picture  leads  to  no  potential  energy  profile,  one  would  have  a  quantum- 
mechanical  dipole  layer,  since  the  electron  density  cannot  drop  from  nd  to  zero 
on  a  length  scale  shorter  than  Xr  Here  XF  is  the  Fermi  wavelength  of  the 
degenerate  electron  gas  inside  the  doped  regions  of  the  structure.  A  “back  of 
the  envelope”  calculation  based  on  this  idea  leads  to  a  quantum  mechanical 
potential  well  in  the  positive  background  region  of  the  structure  whose  depth 
is  comparable  to  the  depth  of  the  thermally  induced  space-charge  well  seen 
by  holes  at  room  temperature.  This  implies  that  a  full  quantum- mechanical 
treatment  might  be  required  in  order  to  obtain  reliable  space-charge  fields 
(and  hole  binding  energies)  in  these  DBRT  and  similar  structures.  A  less 
important  approximation  involved  replacing  the  Fermi  distribution  by  a 
Boltzmann  distribution  to  obtain  equation  (1).  Although  conditions  required 
for  this  approximation’s  validity  are  violated  to  some  extent  in  the  structure 
we  considered,  corrections  due  to  electron  degeneracy  are  small  (in  the 
structure  we  considered)  and  do  not  change  the  hole  binding  energies 
substantially. 

In  view  of  the  unsatisfactory  nature  of  some  approximations  made  in  this 
work,  we  are  presently  applying  density  functional  theory  to  this  problem  and 
will  report  on  that  work  in  the  future. 
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